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FOREWORD 


This  Quarterly  report  was  prepared  by  RCA  Laboratories,  Princeton,  New 
Jersey  under  Contract  No.  F30602-71-C-0356.  It  describes  work  performed  from 
9  December  1971  to  8  March  1972  in  the  Communications  Research  Laboratory, 

Dr.  K.  H.  Powers,  Director.  The  principal  investigator  and  project  engineer 
is  Dr.  D.  A.  de  Wolf. 

The  report  was  submitted  by  the  author  on  7  April  1972.  Submission  of 
this  report  does  not  constitute  Air  Force  approval  of  the  report's  findings  or 
conclusions,  it  is  submitted  only  for  the  exchange  and  stimulation  of  ideas. 

The  Air  Force  Program  Monitor  is  Lt.  Darryl  P.  Greenwood. 
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This  technical  report  has  been  reviewed  and  is  approved. 


ABSTRACT 


Several  statistics  of  the  focal  spot  of  a  focussed  laser  beam  In  turbulent 
air  are  computed.  The  average  area  Is  recomputed  for  horizontal  propagation, 
and  corrected  graphs  are  presented.  Comparison  to  the  half-power  area  Is  made, 
and  It  Is  concluded  that  the  irradlance-weighted  radius  variance  Is  a  slightly 
better  criterion,  although  the  differences  are  small.  The  calculations  are 
extended  to  slant  paths.  The  power  spectrum  of  the  area  scintillation  Is  shown 
to  be  proportional  to  that  of  angle  of  arrival  of  a  ray.  It  Is  computed  for 
frozen  flow  and  for  random  flow  with  zero  mean  velocity.  Both  spectra  are 
Initially  flat,  then  decay  as  f”2/3.  A  composite  formula  is  suggested  for 
the  general  case. 


Ill 


SUMMARY 


Atmospheric  turbulence  causes  the  focal  spot  of  a  focussed  laser  beam  tc 
scintillate,  i.e.,  the  irradiance  distribution  in  the  focal  plane  varies  irreg¬ 
ularly.  In  the  first  quarter  of  the  contract  period  (9  June  to  8  September  1971), 
a  simple  calculation  of  the  variance  of  an  irradiance-weighted  radius  of  the 
focal  spot  was  presented.  It  yielded  an  average  focal-spot  area  equal  to2the 
sum  of  the  diffraction-limited  area  and  an  area  proportional  to  L  <m  '  Cn 
(L  -  range,  ■  microscale  wavenumber,  Cn2  *  refractive-index  structure 
constant). 

There  exist  arguments  whether  such  a  variance  represents  an  effective  area 
as  well  as  the  half-power  area  in  the  focal  plane.  Moreover,  some  numerical 
errors  crept  into  the  analysis  and  the  graphs.  Hence,  the  entire  matter  is 
presented  again  in  Section  I  with  errors  rectified,  and  with  a  considerable 
discussion  of  the  "differing"  definitions  of  effective  focal-spot  area.  We 
conclude  that  the  differences  are  unimportant,  and  that  the  irradiance-weighted 
variance  is  more  proper. 

In  Sections  II  and  III,  the  analysis  is  extended  to  slant  paths,  ground-to- 
air,  and  air-to-ground.  The  theory  summarized  by  Wvngaard,  et  al.[5]  is  uti¬ 
lized  for  incorporating  the  height  dependence  of  Cn  (z)  into  the  theory.  The 
effects  are  much  less  than  in  ground- level  horizontal  propagation,  and  weakest 
for  the  ground-to-air  case.  Graphs  are  presented  that  allow  one  to  convert  the 
horizontal-propagation  parameters  into  slant-path  equivalents,  when  altitude  of 
laser  or  receiver  is  known. 

In  Section  IV,  the  power  spectrum  of  focal-spot  radius  scintillation  is 
computed.  It  is  proportional  to  the  power  spectrum  of  angular  fluctuations  of 
a  ray.  We  compute  the  latter  for  two  cases:  (i)  Taylor's  hypothesis,  i.e., 
frozen  flow  with  a  steady  wind  velocity,  and  (ii)  random  velocity  distribution 
with  zero  mean.  An  approximate  formula  is  constructed  for  the  composite  of 
these  two  cases:  it  yields  a  spectrum  that  is  a  function  of  a  square  velocity 
v2  >  Ut2  +  16AU2/3.  Here,  Ut  is  the  mean-velocity  component  across  the  beam, 
and  AU*  is  the  velocity  variance.  The  spectrum  is  flat  when  w  <  v/L0  (w  is  the 
angular  frequency  variable  of  the  spectrum,  L0  is  the  macroscale) .  It  decays 
as  u_2/3  for  v/L0  <  w  <  Kmv,  and  cuts  off  quite  sharply  when  w  exceeds  umv. 
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I.  HORIZONTAL  PROPAGATION  PATH 


In  the  First  Quartetly  Report[l],  a  formula  was  derived  for  the  effective 
focal-spot  radius  tl  of  a  focussed  laser  beam  propagating  through  turbulent  air 
characterized  by  strength  Cn2  and  scales  lQ  *>  5.92/icm  and  L0.  The  definition 
of  is 


rL2  *fd2 p  p2  <I(?,L)>//d2p  <I(p,L)> 


(1) 


hence  it  is  an  average-irradiance  weighted  variance  of  the  distance  in  the  focal 
plane  from  thje  central  axis.  The  result  we  obtained  is 


r  2  =  (L/kr  )2  +  1.3  L3k  1/3C  2 

Li  O  — 


m 


(2) 


for  a  Gaussian  beam.  The  general  result  is 


r  2  =  r  2  +  1.3  iA  1/3C  2  ' 
L  Lo  m  n 

i  '  : 


(3) 


where  rLo  is  defined  by  (1)  if  <  I(p,l1)>  is  replaced  by  the  free-space  I0(p»L) 
The  result  (2)  is  essentially  the  same  as  that  of  Varvatsis  and  Sancer[2]. 

1  :  . 

Lutomirski  and  Yura[3]  define  another  effective  focal-spot  radius  tl' 

(z0  in> their  notation)  which  is  defined  as  a  half-power  radius  such  that 
<iI(^’,L)>  ■  1/2  <  J(0,L)>  for  p  ■  tl’.  They  obtain 


r  '  *  0.9  L3/2  k  1/6C 
"  !  m  n 


for  0.4k2LJl  5/3C  2  »  1 
'  s  o  n 


*  0.7L8/5k1/5C  6/5  for  0.4k2U  5/3C  2  «  1 
i  n  on" 


(4) 


There  is  aiso  a  lower  boupd  on  L  for  the  second  criterion  in  (4)  but  it  is  of 
the  order  of  one  or  sleveral  mSters  for  systems  with  apertures  of  at  least  10  cm 
diameter,  and  hence  not  of  much  practical  Importance. 

!  :  : 

i  The  results  (3)  and  (4)  appear  different  for  L  less  than  (0.4k^£o^^Cn2)”'*', 
a  distance  of  roughly  7  km  at  0.6  pm  and  more  than  2000  km  at  10.6  pm.  This 
difference  thus  appears  to  manifest  Itself  in  practical  situations,  particularly 
in  the  infrared  regime.  At  the  tjLme  of  writing  -  early  1972  -  there  appears 
to  be  aisagreement  among  workers  in  the  field  about  relevance  and  applicability 
of  results  (3)  and  (4).  The  purpose  of  this  section  is  to  clarify  the  matter. 

I  i 

»  !  * 

The  peculiarity  of  result  (3)  is  the  following:  for  truncated  beams  - 

particularly  for  the  so-called  Airy-disk  case  -  one  finds  rT  ■  hence  rT  ■  <*> 

Lo  L 

:  I  I 


1 


I 


depends  on  the  micro- 


and  the  result  appears  meaningless, 
scale  ZQ  (or  tc  “*) . 


Moreover , 


2 


The  peculiarity  of  result  (4)  lies  in  the  somewhat  arbitrary  definition 
of  as  a  half-power  radius,*  and  the  dependence  of  rj/  upon  the  macroscale  L0. 


We  wish  to  make  a  case  for  the  statements  that  result  (3)  is  more  funda¬ 
mental  than  result  (2)  -  also  in  the  case  that  rL0  -*■<*>-  and  that  the  differ¬ 
ence  between  (ri^-rLo  )^'2  and  rj/  is  rather  small  for  all  practical  cases. 


In  °rder  to  do  so,  it  is  first  noted  that  the  definitions  of 
(rL  -rLo  2  =  rLB*  and  1*25  rj/  coincide  when  0.4k2L£c>  ^Cn2  >>  1.  The 
reason  is  that  the  irradiance  profile  at  z  *  L  has  become  Gaussian  in  which 
case  the  half-power  radius  differs  from  the  irradiance-weighted  variance  by 
only  a  numerical  factor  ^  1.25. 


The  next  point  is  that  (3)  is  a  very  general  result,  as  noted  on  page  9 
of  Reference[l] .  In  Appendix  I,  we  give  a  derivation  well-known  in  radio 
astronomy [4]  where  an  analogous  problem  exists  in  analyzing  the  power  spectrum 
of  backscattered  radio  pulses.  The  spectrum  is  a  convolution  of  the  spectra  of 
the  target  scattering  amplitude  and  the  pulse  envelope.  Hence  the  variance 
width  of  the  spectrum  (which  tells  you  something  about  the  Doppler  content  of 
the  astronomical  object  viewed)  is  the  Pythagorean  sum  of  the  variance  widths 
of  the  two  component  spectra  in  the  convolution.  The  variance  width  of  a  square- 
pulse  spectrum  is  infinite.  Hence  the  analogy.  The  difficulty  is  circumvented 
by  regarding  the  difference  of  the  variances  which  can  be  computed  in  any  prac¬ 
tical  case  by  exploiting  the  finite  bandwidth  of  the  observed  spectrum.  In 
our  case,  that  means  that  rjjj2  -  rj^rn^  can  be  computed  in  practice  because 
rLc>  is  not  really  infinite:  it  is  limited  by  finite  signal- to-noise  ratios  and 
if  these  are  taken  into  account  consistently  for  ^  and  rL0  there  is  no  real 
difficulty.  I.e.,  one  computes  from  the  data. 


/  dZp  P2  [<I(p,L)>  -Iq 
0 


<  I  (p  ,L)> 


(5) 


where  pn  is  a  value  of  the  radius  beyond  which  the  patterns  vanish  into  noise. 

If  the  broadening  due  to  turbulence  occurs  well  within  pn,  i.e.,  if  the  result 

of  the  above  calculation  yields  tlb  <<:  Pn»  then  it  follows  that  the  denominator 

of  (5)  is  adequate  since  it  is  essentially  equal  to  the  total  power  emitted  by 

the  laser.  It  is  quite  true  that  the  second  p-moments  of  <  I(p,L)>  and  I0(?,L) 
are  proportional  to  pn,  but  their  difference  is  not.  The  quantity  r^j2  there¬ 
fore  can  be  calculated  and  has  meaning  even  when  the  variance  width  r^02  of  the 
vacuum  irradiance  pattern  is  infinite.  The  irradiance  pattern  is  thus  broad¬ 
ened  in  turbulent  air  by  convolution  with  the  Fourier  transform  13(9)  “  see 
Appendix  I  -  of  the  mutual- coherence  function  of  a  spherical  wave  field  in 
turbulent  air  divided  by  the  equivalent  vacuum  field.  The  quantity  is 

simply  the  variance  width  of  this  Fourier-transform  pattern  IoGi),  where 


*However ,  definition  (II)  of  Ref.  1  yields  essentially  the  same  result  (private 
communication,  D.  P.  Greenwood). 
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In  my  opinion  this  settles  the  question  of  whether  the  resulting  broadening 
should  depend  upon  the  macroscale  L0  or  ujson  the  microscale  £0.  The  broadening 
is  clearly  determined  by  the  shape  of  Ig(q)  •  The  variance  width  ruj^  of  Ig(^) 
is  finite  and  therefore  clearly  a  valid  measure  of  the  width.  It  is  given  by 
(A-10) ,  hence  it  is  a  function  of  the  microscale  .  It  is  fairly  easy  to  under¬ 
stand  why  the  broadening  depends  more  essentially  upon  the  microscale.  The 
definition  of  Ig(^)  in  (A-3)  makes  it  perfectly  clear  that  the  radius  parameter 
<f  =  (k/L){?  of  the  Fourier-transform  pattern  in  the  focal,  plane  is  the  Fourier- 
conjugate  variable  of  the  difference  coordinate  Ap  in  <  B(1)B*(2)>  Consequent¬ 
ly  the  behavior  of  Ig(q)  for  large  q  is  determined  by  the  small  Ap-behavior  of 
<  B(1)B*(2)>  .  These  are  both  therefore  determined  by  the  microscale  lQ.  On 
the  other  hand,  the  macroscale  L0  determines  the  behavior  of  <  B(1)B*(2)>  for 
large  Ap,  hence  of  Ib(<J)  for  small  q. 

What  we  are  saying  is  in  effect  that  the  half-power  radius  ri/  is  not 
really  a  good  measure  of  the  broadening  (although  we  will  show  further  on  that 
it  isn't  bad  either).  Remember  that_^the  pattern  <I(p,L)>  in  the  focal  plane 
is  a  convolution  of  IQ(p,L)  with  Ib(p,L)  =  IgC$).  When  In  is  broad  compared 
to  I0,  obviously  the  broadening  is  given  by  rLB^  ^  L\m^/3cn^  and  the  first 
criterion  for  rj/  given  in  (4)  must  hold.  When  IB  is  narrow  compared  to  IQ 
the  broadening  is  also  given  by  tlb  ^  L^Km^/3(jn2f  but  this  quantity  is  then 
small  compared  to  the  width  of  the  main  lobe  of  I0,  and  therefore  not  interest¬ 
ing.  The  interesting  case  is  when  ruj  is  comparable  to  the  main- lobe  width  in 
which  case  Lutomirski  and  Yura[3]  utilize  ri/  ^  18/5^1/5^6/5  an(j  the  second 
criterion  in  (4).  At  infrared  frequencies,  most  propagation  paths  of  interest 
fall  into  this  case  when  Cn2  lO-!^-^'-*.  what  happens  in  this  case  is  that 
the  nulls  between  main  lobe  and  side  lobes  of  I0  are  rapidly  filled.  The  total 
patterns  <  I>  given  in  [3]  indicate  this  at  p'ths  L  larger  than  several  hundred 
meters  (L  -  100  zc) .  The  half-power  radius  then  becomes  rather  arbitrary. 

There  is  another  way  to  examine  this  difference  between  r^/  and  rjjj.  It 
is  simply  to  compare  them!  In  fact,  let  us  compare  tlb  and  1.25  ri/  to  each 
other,  since  these  coincide  for  a  Gaussian  beam. 


1.25  r’  L8'V/5Cn6/5 

^  L3/2 ;  1/— 
m  n 


LB 


(■ 


<  C 

m  n  / 


2\l/10 


(6) 


For  fixed  frequency,  this  ratio  varies  with  (LCn  )  ,  i.e.,  very  weakly. 

Clearly,  the  definitions  are  not  very  different.  More  serious,  however,  is 
the  restriction  of  validity  of  the  numerator  of  (6)  to  the  condition  (p.  1657 
of [3]), 


1  2 

10  <  0.4k  L  L 


5/3 


n 


<  10 


(7) 


for  errors  of  less  than  25%.  By  taking  the  10-th  root  of  (7),  we  obtain 


3 


2 

A  typical  large  surface  value  for  tcmL0  is  10  ,  hence  we  note  th-it  in  the  region 
of  validity  of  the  expression  for  r  ' ,  we  have 

Lt 


0.4  <  1.25rL  /rLB  <  0.9  (9) 

Relationships  (6)  and  (9)  indicate  clearly  that  there  is  not  much  numerical 
difference  between  r^’  and  r^g  in  the  region  of  validity  of  the  cited  expression 
for  ri/  .  On  the  other  hand,  the  expression  for  r^g  is  valid  over  a  wider  range. 
Hence  it  is  concluded  that  (3)  is  to  be  preferred  to  (4)  ,  although  careful  use 
of  (4)  does  not  lead  to  appreciable  differences . 

2 

Horizontal  propagation  is  characterized  by  constant  Cn  along  the  propaga¬ 
tion  path.  It  is  therefore  fairly  simple  to  plot  (3)  for  this  case.  Figures 
1  and  2  of  Reference[l]  contain  some  inaccuracies,  hence  new  graphs  are  sup¬ 
plied  in  this  report. 


Figure  1. 

The  free-space  focal-spot  radius 
rLo»  and  its  average  broadening 

rLB  «rL2>  =  rLo2+rLB2)  as  Unc¬ 
tions  of  range  L  for  lO-1**  <  Cn^ 
in  m“2/3  <  10“12. 
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In  Figure  1,  we  plot  rj,  vs  L  according  to  (3)  .  The  solid  lines  represent 
rLB»  the  second  term  in  the  right-hand  side  of  (3).  The  dashed  lines  represent 
r^o  for  two  different  frequencies. 

In  Figures  2  and  3,  we  plot  roc  vs  L.  The  critical  laser  aperture  roc  is 
defined  to  be  such  that  r^o  =  If  roc  *8  increased,  then  r]_0  is  decreased 

and  rj,  %  r^g.  That  indicates  that  larger  apertures  than  those  with  radius  roc 
will  not  yield  smaller  average  focal  areas:  a  lower  limit  irrLB^  masks  the 
diffraction  area.  By  setting  r^o  =  L/kroc  equal  to  rig  we  obtain 


O.l  I  10 

L  in  km  - p- 


Figure  2.  Maximal  aperture  radius  roc  required  to  obtain  a  minimal 
average  focal-spot  area  as  a  function  of  range  L  for 
10~14  <  CnZ  ln  m-2/3  <:  io~12  at  X  -  0.6  pm. 
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Figure  3.  Maximal  aperture  radius  roc  required  to  obtain  a  minimal 
average  focal-spot  area  as  a  functlqn  of  range  L  for 
10'14  <  cn2  in  m'2/3  <  i0-12  at  A  _T10.6  um> 
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II.  SLANT-PATH  PROPAGATION :  GROUND-Td-AIR, 


In  this  section,  We  will  modify  (^)  for  the  case  of  slant-path  propaga¬ 
tion  in  which  case  Cn,2  is  a  function  of  altitude.  To  do  so,  we  note  that  the 
main  difference  with  horizontal  propagation  is  tha,t  instead  of  obtaining  in 
(A-10)  a  factor 


00 

/ 


Lc"  /  dKK  4>(K)  ;  ' 

i 

P 

from  the  limit  p  +  0  in  (A-9) ,  we  should  be  obtaining 


lim 

p-*0 


.  -2 
4p  i 


L  r 

f  dse2(s)  J 


dKK4> (K)  [  1  -  J  (Kps/L)] 


(11) 


1 

/ 


,  dse2(s)  (s/L)2  1  j  dKK3$(K) 


for  a  slant  path  characterized  by  a  path  parameter  s.  Equation'  (11)  inserted 
ilnto  (A-10)  yields 


rj-  3,9  A 


3/3  f  ds(s/L)2  C2(z) i 
m  j 'i  n 


(12) 


In  adapting  (.12)  to  a  situation  where  the  laser  is  located  at  an  altitude  zQ 
close  to  the  ground  and  pointed  upwards  at  an  elevation  'angld  of  0  with  the 
horizop  towards  an  objeqt  at  distance  L,  it  is  noted  that  path  coordinate,  s  is 
related  to  altitude  z  through  s  ■  (z-z0) /sin0. 

i  ■ 

2 

The  behavior  of  Cn  (z)  with  altitude  is  fairly  complicated  -  particularly 
at  altitudes  above  1  km.  The  Hufnagel  curve,  Figure  4,  indicates  the  diffi¬ 
culties:  i occurrence  of  random  disturbed  layers  and  dependence  of  lower- 
altitude  Cn  upon  time  of  fday.  Wyngaard ,  et  al.  have  surveyed  recent  develop¬ 
ments  in  altitude  dependence  of  Cn2[5] .  Their  results  yield. 


Cn  <2)l-  Cn  (*0)(VT)2(2o/2)i/3i<1+T/‘s>T2/3 


x  exp[-2(z-zQ)/h]  , 


(13) 


«  *7* 


i  ■ 


i 


i 


rt 

s. 


figure  4.  Cn(z)  in  cm  vs  altitude  z;  the  "Hufnagel  curve” • 
The  Wyngaard,  et  al.  predictions  [Reference  5]  are 
included . 


where  T  is  the  temperature  at  altitude  z  (T0  at  zG) ,  h  is  the  scale-height  of 
pressure  at  sea  level  (h  «  8.3  km),  and  fcs  is  a  length  parameter  related  to  the 
Monin-Obukhov  length  L^g  by  is  ■  -Ljjo/7.  This  model  predicts  (for  2z  «  h) . 


r  2  a  ,-4/3 

C  a  Z 

n 

sunny  day  -  midday  (£ 

c  2  .  *-2/3 

n 

dawn,  dusk  (£  -+  °°) , 

s 

and  Figure  4  bears  the  prediction  out  quite  well.  The  older  Fried  model, 
cited  by  Titterton[6] ,  and  Sutton[7]  will  therefore  not  be  followed.  Note 
that  Equation  (12)  is  a  function  only  of  the  receiver  altitude  Z  «  LsinG,  and 
that  it  can  be  rewritten  as 
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(14) 


4  <z>  -  4  <°>-F<z> 

z 

F(Z)  -  |  j  dz  exp(-2z/h)(z/Z)2 
z 

o 

x  (1+z/z  )  2/3[l+(z+z  )/£  ]  2^3 

O  U  D 

Note  that  the  temperature  dependence  on  altitude  has  been  ignored.  It  is 
variable,  but  weak  compared  to  the  other  s-dependent  parameters  in  the  F(Z) 
factor  of  (14) .  Certain  approximations  simplify  the  calculation  of  (14)  : 

(i)  sunny  day  -  midday  situation :  We  may  replace  the  last  two  fac¬ 
tors  in  (14)  by  (z0*s)2/3  (z)"4'3  provided  Z  »  z0,£s.  For  Z  «  zD 
we  obtain  F(Z)  *=  1.  The  critical  altitude  Zc  for  which  Zc  -  zD,  !LS 
is  of  the  order  of  1  m  or  less  when  z0  88  1  m,  and  L  >  1  km.  Hence, 
slant  paths  are  handled  easily  for  most  non-zero  angles  0  by 
means  of  this  simplification. 

(ii)  dm.m}  dusk  situation".  Because  l8  °°  in  this  case,  we  replace 
the  last  two  factors  in  (14)  by  (z/z0)~2/3  under  the  proviso 
that  Z  »  zQ,  which  is  always  true  for  any  airborne  target 
under  consideration. 

The  remainder  of  this  section  is  devoted  to  straightforward  numerical 
expression  of  (14) .  It  is  easily  seen  from  (10)  and  (14)  that 

roc<Z)/roc(0>  -  rLB(0)/rLB<Z)  •  tF<Z)]-1/2.  (15) 


Hence,  it  suffices  to  plot  (15)  as  a  function  of  Z  in  order  to  infer  the  slant- 
path  quantities  from  those  at  zero  altitude  for  an  equivalent  path  of  length  u. 


(i)  sunny  day  -  midday".  With  the  previously  discussed  simplifica¬ 
tions,  it  follows  that 

F(Z)  *  1.8  [(z  i.  )1/2/Z]4/3  for  Z  «  h/2 
o  s 


*  0.85  h5/3(z  i  )2/3/Z3  for  Z  »  h/2  , 
o  s 

the  intermediate  regime  is  governed  by  the  incomplete  gamma 
function  from  which  (16)  is  derived. 

(ii)  dam t  dusk:  Likewise,  we  obtain 


(16) 
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for  Z  «  h/2 
for  Z  »  h/2. 


(17) 


F(Z)  *  1.3(z  / Z) 
o 


2/3 

2/3/2.3 

o 


In  our  view,  the  Z  >>  h/2  forms  in  (16)  and  (17)  should  not  be  regarded 
very  seriously  because  they  depend  upon  validity  of  the  model  (13)  above 
2  =  h/2  %  4  km,  and  that  is  certainly  above  the  altitudes  where  this  model 
is  valid. 

In  Figure  5,  we  plot  the  Z  «  h/2  forms  of  [F(Z)]"1^2  vs  Z,  using 
*  z0  ss 1  m.  The  plots  are  substantially  the  same  as  those  in  the  Second 
Quarterly  Report  (9  September  -  8  December  1971,  sequel  to  Reference [1]) . 
However,  it  is  simpler  to  plot  the  ratios  as  functions  of  altitude  Z,  and  to 
modify  the  horizontal— path  results  of  Figures  1  to  3  by  these  altitude  factors 
when  the  image  at  the  end  of  path  L  is  at  altitude  Z. 


Figure  5.  Modification  factors  for  tl0  and  roc  (from  Figures  1  to  3) 
for  ground-to-air  propagation  for  laser  transmitters  at 
altitude  Z. 
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Clearly,  the  effect  of  turbulence  Is  much  less  pronounced  on  slant  paths. 
It  should  be  stressed  that  the  curves  In  Figure  5  depend  on  and  on 

Is1/3,  and  that  we  have  chosen  z0  **  1  m,  £s  w  1  m.  The  slant-path  results  of 
Figure  5  should  be  scaled  accordingly,  If  different  values  for  zQ  (laser 
height)  and  £s(*  -Lmq/7)  are  utilized.  Note  that  Z  ■  L  sin  6;  these  curves  are 
functions  only  of  the  product  L  sin  6.  Each  point  describes  a  variety  of  path- 
lengths  L  and  elevation  angles  6  subject  to  L  sin  0  ■  Z. 
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III.  AIR-TO-GROUND  PROPAGATION 


If  it.  is  the  transmitter  that  is  at  an  altitude  z  ■  Z,  and  the  receiver 
at  z  -  zQ  such  that  the  propagation  path  has  a  length  L  at  an  inclination  of 
-0  with  the  horizon  (ignoring  curvature  effects),  then  a  slightly  different 
calculation  ensues.  Instead  of  (12),  we  obtain 

4  ■  3‘9  /  d8<1  -  S/L)2  °n(2)  ' 

0 


Because  the  dominant  part  of  the  integral  is  in  the  first  150  m  or  so  above 
the  ground,  where  (1-s/L)2  *  1  for  altitudes  Z  of  500  m  or  more,  and  because 
the  model  for  Cn2(z)  is  not  trustworthy  above  several  hundred  meters,  we  may 
just  as  well  ignore  the  factor  (1-s/L)2  in  (16)  and  adapt  (14)  to  the  present 
case  by  writing 

r?„(Z)  =  r2  (0)  F*  (Z) 


F*  (Z) 


L 

if 


(19) 


dz  exp(-2z/h)  (1+z/zq)  2^[l+(z+zo>/&g]  2^3 


We  have  utilized  Lsin0  *  Z  in  this  formula,  as  in  the  previous  section.  The 
numerical  work  is  a  little  more  complex  than  in  Section  II. 


(i)  Dawn t  dusk:  In  this  case  (17)  simplifies  to 


F'(Z) 


Z/zo 

J  dx  (l+x)"2/3  exp(-2xZ/h), 
1 


(20) 


and  when  Z  «  h/2  (we  need  not  consider  the  other  case  because 
the  model  is  then  too  uncertain)  it  follows  that,  to  good 
app  r oximat ion , 


F'(Z)  *  9(z  /Z)2/3  for  Z  «  h/2. 

(ii)  Sunny  day-midday:  To  sufficiently  accurate  approximation,  (19) 
can  be  simplified  to 


F'(Z) * 


3zl/3a2/3  ZHo 
o  s 


j  dx  x"2^3(l+x)"2^3exp(-2xzo/h) , 


(21) 


(22) 


and  a  further  approximation  is  to  replace  Z/zQ  in  the  upper 
bound  of  integration  bv  ».  Replacing  x  by  y-1/2,  we  obtain 
an  integrand  (y2-l)-2/2,  in  which  we  may  ignore  the  -1  term; 

i  *6  •  i 
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F'(Z)  ^ 


(23) 


.  1/3.  2/3  2 

3--°-  - - /  dy  y'4/3  *  9  *  1/3l  2/3/Z. 

ZJ  os 

3/2 

The  resulting  expressions  (21)  and  (23)  are  Inserted  Into  (15)  which  Is 
plotted  In  Figure  6  for  air-to-ground  propagation.  The  effects  of  turbulence 
are  somewhat  larger  In  this  case  than  in  the  equivalent  opposite  case:  ground- 
to-air.  The  reason  Is  that  the  filter  factor  of  (12)  eliminates  contributions 
from  low-altitude  Cn2(z),  whereas  that  of  (16)  removes  hlgh-altltude  (much 
lower)  Cn2(z).  This  may  appear  surprising,  but  we  must  remember  that  the 
focussed  bundle  consists  of  spherical  wavelets  leaving  the  aperture  at  z  *  0 
and  coming  together  close  to  z  *  L,  each  at  its  own  angle.  Turbulence  affects 
the  direction  of  a  spherical  wave  most  when  the  distance  from  the  point  of 
transmission  Is  large.  The  effect  Is  analogous  to  the  "dirty  windshield" 
effect:  a  driver  with  a  dirty  windshield  can  be  blinded  by  the  headlight  of  an 
oncoming  car,  but  the  oncomer  cannot  see  that  the  windshield  of  the  other  car 
is  dirty. 


r0C(Z)  /roc(O)  ■  rLB(0)/rLB(Z) — ► 

Figure  6.  Modification  factors  for  rL0  and  roc  (from  Figures  1  to  3) 
for  air-to-ground  propagation  for  laser  transmitters  at 
altitude  Z. 
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IV.  SCINTILLATION  RATES:  ANGULAR-FLUCTUATION  POWER  SPECTRUM 

l 

A.  INTRODUCTION 

i 

t 

So  far,  we  have  discussed  only  the  average  focal-spot  broadening  due  to 
atmospheric  turbulence.  What  really  happens  is  that  the  focal  spot  scintillates 
around  the  average.  Consequently  the  next  questions  deal  with  higher-order 
statistics  such  as  the  variance  of  the  focal-spot  area  around  the  mean,  and 
with  the  time  scales  on  which  this  phenomenon  occurs. 

i 

We  will  occupy  ourselves  in  this  section  with  the  time  scale  of  scintilla¬ 
tion;  i.e.,  with  scintillation  rates.  There  are  many  ways  to  define  a  signifi¬ 
cant  measure  of  scintillation  rates,  but  the  one  most  cra™ionly  used  is  the 
power  spectrum  of  a  scintillating  quantity.  The  power  spectrum  is  the  Fourier 
transform  (with  respect  to  separation  time  x)  of  the  autocovariance  of  the 
scintillating  quantity  with  itself  at  a  time  x  later.  It  is  a  function  of  the 
conjugate  variable  u,  an  angular  frequency.  i  1  !.  i 


Let  us  try  to  define  an  appropriate  au  toe  ova  ri  ancle  function  Rtr(x)  that- 
reduces  to  rLB^  when  x  -*■  0.  To  do  so,  we  note  that  the  average  of  Equation 
(A-l)  is  a  special  case  of 

< E(r;t)E*(r;t+x)>  i 

/_  1  (2 
*\fr292  <B(r,rlit)B*(?,?2it+T)>Uo(r1)U*(;2)exp[lk$.(?1-J2)/L]. 

i 

When  x  -*■  0,  we  note  that  the  above  formula  reduces  to  tha(:  for<I(r)  >.  Soi 
the  new  definition  that  we  introduce  for  RLB(x)  is  obtained  from  (5;  by  re¬ 
placing  <1  (p,L)  >  in  the  numerator  by  (24),  viz., 

oo  m 

RLB(t)  "  /  d2pp2[  <E(?;t)E*(r;t-K)>  -Io(r)]/  f  d2p  <!(?)>  ,  J  (2 


where  r  ~  (p ,L) .  In  f act,  the  numerator  of  this  definition  can  be  obtained  by 
multiplying  p  [E(r;t)-EQ(r;t) ]  by  the  complex  conjugate  of  p  [E(r;t+x)-EQ(r;t+f) ] , 
averaging  this  and  then  integrating  over  the  focal  plane  (the  cross  terms  are 
zero  on  the  average) .  It  is  therefore  plausible  that1  R(x)  is  physically, 
analogous  to  the  auto covariance  of  a  focal-spot  radius  at  time  separation  x. 

Let  us  therefore  adopt  (25)  as  the  desired  autocovariance.  i 

....  th®  derivation  in  Appendix  I  is  followed  with, this  one  difference:  i 

t0  be  used  instead  of  (A-2) ,  it  Is  readily  seen  that  we  obtain,  instead 
or  (A-10) ,  !  1 


i 


i 


J 


I 


*lb<*>  ■  5T  /  dKK3  *41  (KiT)  • 

:  0  , 

i  •  : 

O  I  !  1  4.  4. 

where  e  $^(K;t)  is  the  Fourier  transform  of  <  6e (ri;t)6e (rojt+x)  > 
also  over  t,  with  respect  to  thje  difference  coordinate  ri-r^  ■  Ar. 
noted  by  comparing  (26)  with  Appendix  II  of  Reference  [1]  that 


(26) 


averaged 
It  will  be 


Rlb (t)  “  J"  <  6^(t)-6^(t+x)>  , 


(27) 


I  ;■ 

i.e.,  that  thq  autocovariance  of,  a  weighted  focal-spot  radius  is  intimately 
associated  with  the  autocorrelation  of  the  angular  deviation  6$  of  a  ray  over 
a  path  of  length  L.  For  those  who  feel  uncomfortable  with  the  definition 
given  for  RlbCt)»  we  can  therefore  state  that  another  physical  interpretation 
of  this  quantity  is  that  it  is  proportional  to  the  autocovariance  of  the 
angular  deviation  of  a  ray.  Let  us  therefore  define 


1  VT) 


fiW  ♦*«•*>•  ! 


(28) 


■  1  1 

The  qonnectlon  between  Rg  and  Rlb  follows  from  (27) .  We  shall  now  restrict 
ohrselves  to,  (28) ,  and  to  horizontal  propagation  with  uniform  Cn2  along  the 
path  for  simplicity.  Extensions  to  slant  paths  will  be  pointed  out  where 
germane.  .  ’  ' 


The  major  point  in  further  developing  (28)  is  to  make  a  good  hypothesis 
for  $4(K;t).  By  definition  , 


M*;t) 


■/ 


d3Ar 


<5e(r1;t) 


6e*(r2;t+x)>  e1K  Ar, 


(29) 


where  averaging  oyer  t;ime  is  implied.  In  Appendix  II,  we  utilize  Favre's 
hypotheses [8]  to  find  (A-17)  as  the  major  step  in  ensuihg  development.  In¬ 
serting  (A-17)  into  (28),  one  obtains,' 

00 

i  2  j  /•  1  .  . 

Re(T)  -|J-  J  dKK3  $(K)  Jo(KUt  t)  exp(-  y  K2AU2t2) ,  (30) 

1  ,  !  0 

where  Up  is  the  component  of  the  mean  velocity  normal  to  the  propagation 
direction,  and  AU2  is  the  variance  of  the  wind  velocity  In  isotropic  turbulence. 
When  AU2  ■  0,  Equation  (28)  describes  "frozen  flow"  or  "Taylor’s  hypothesis": 
the  motiion  is  a  steady  convection  'of  a  random'  spatial  structure  of  eddies 
(that  are  immobile  with  respect  £o  each  other)  across  the  beam  with  velocity 


15 


1 


component  Uf  By  using  Equation  1-3  of  Reference  [1],  and  transforming  K  to 
a  new  variable  x  ■  r2l0^,  we  rewrite  (30)  as 


V’> 


15 


16ir 


^  f  dxx(l+x)"11/6  Jo(o>tt 


x  exp[-x(A(D2T2/4  +  1/k2  L2)  ] , 


(31) 


where  uj.  *  Ux/Lq  and  Aw2  -  16AU2/3LQ2  are  the  frequencies  obtained  by  dividing 
the  previously  introduced  velocities  by  the  macroscale  LQ.  Finally,  the 
potter  spectrum  of  angular  fluctuations,  W0(w),  is  the  Fourier  transform  of  (31): 

00 

W0(«)  =  2  /  dt  R0(t)  cos  wt  (32) 

0 

A  simplification  can  be  made  by  noting  that  x(l+x)  =  (1+x)  ^^-(1+x) 

It  allows  one  to  rewrite  (31)  and  (32)  as 


WQ(u>) 


15.7 

8ir 


2. 

[W0 (ai.5/6)  -W0 (oi.ll/6)]. 


(33) 


with  the  following  definition  of  W0(oi,q): 

00  00 

WQ (oi,q)  •  J"  dx  J  dx  (1+x)  q  J0(<*>tt  *^c) 

0  0 

x  exp[-x(Aoi2T2/4  +  1  /k2  L2)]  cos  oit  . 


(34) 


Thus,  we  attempt  to  compute  (34)  for  q  ■  5/6  and  for  q  *  11/6,  and  insert  the 
difference  into  (33)  to  obtain  the  power  spectrum.  Unfortunately,  the  double 
integral  in  (34)  does  not  appear  integrable  into  a  tabulated  function  in 
general.  However,  we  can  do  two  special  cases,  namely  Aw  ■  0  and  ux  ■  0. 

These  two  cases  are  the  subject  of  subsections  B  and  C. 

B.  FROZEN  FLOW  (TAYLOR'S  HYPOTHESIS),  AU2  -  0 

In  this  case,  it  is  easiest  to  carry  out  the  Fourier  transformation  in 
(34)  first,  inverting  the  order  of  integration,  to  obtain 

00 

-  WQ  («,q)  «  y  dx  (l+x)-q(xui2-ui2)  1^2  exp(-x/K2  L2) .  (35) 

2,  2 

0)  /wT 
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Let  us  also  Introduce  the  frequency  Of  ■  KjnU^;  it  corresponds  to  a  velocity 
divided  by  a  microscale.  Equation  (35)  is  easily  transformed  into, 


WQ(u>,q) 


U)( 


2q-2 


r“2/nT 


w 

/ 


.  -1/2.  .  2 

dyy  (y+ui 


2.-q 

V 


-y/K 


(36a) 


The  second  form  was  obtained  by  utilizing  Equation  (14.4.12)  of  Reference  [9], 
and  converting  the  Whittaker  function  in  the  result  to  a  Kummer  function 
U(a,b,z)  as  in  Equation  (13.1.33)  of  Reference  [10].  Let  us  see  what  W0(w,q) 
looks  like  for  different  parameter  regimes : 

(i)  w  «  ftj.: 

For  the  case  that  to  is  much  less  than  the  microscale-defined 
frequency  U«p/i0»  we  can  utilize  Equations  (13.5.10)  and 

(13.5.12)  of  Reference  [10]  to  obtain 


W0(a>,q)  * 


r(q-i/.2_) 

r(q) 


1/2-q 


(36b) 


(ii)  (o  »  flT: 

In  this  case,  we  utilize  the  asymptotic  form  of  Rummer's  func¬ 
tion  U(a,b,z)  -*■  z_a  to  find, 


W0(<o,q) 


(36c) 


This  formula  is  not  very  interesting  since  u>  »  0-p  implies  a  strong 
Gaussian  decay  factor  making  (36c)  quite  negligible  compared  to  (36b). 

Note  that  the  power  spectrum  is  flat  below  b>  «  Uj/L0,  and  then  decays  as 
w“2/3  for  Uj/x0  <  (o  «  Ul/1^.  As  u  - approaches  Vj/i0,  the  decay  sharpens -into 
a  Gaussian  cut-off,  but  this  has  as  much  meaning  as  the  Gaussian  cut-off  in 
*(K). 
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C.  RANDOM  VELOCITY  (NO  MEAN),  UT  -  0 

Returning  to  (34),  we  again  carry  out  the  Fourier  transformation  first 
to  obtain. 


We("’q)  =  ^  /  d*x"1/2a+xrq  exp  j  )  *  <37> 

«  \xAw  ic  L  / 

0  '  m  of 

This  case  is  slightly  harder  to  handle  than  the  previous  one.  We  can 
drop  the  x/(cm2L02  term  in  the  exponential  because  the  first  term  in  it  will 
take  care  of  convergence.  With  tba  simplification,  we  set  x  -  1/y  in  (37)  to 
obtain 


WQ(w,q)*>  ^  J  dyyq“3^  (l*y)“q  exp(-yw2/Aw2) 
0 


F (q~l/2) 


—  u(q  - 

Aw  y  2* 


(38a) 


The  development  of  the  second  form  is  the  same  as  that  in  (36a) .  It  can  be 
seen  that  the  general  shape  is  similar  to  (36b)  by  using  the  asymptotic  forms 
of  the  Rummer  function  in  (38a) : 


W  (w,q)  'v  -1/2)  JL.  for  u  «  ^w, 

6  r(q)  4“ 


-v  r  (q-1/2) 


/tT  /  w_\ 
Aw  \Aw/ 


l**2q 


for  w  >>  w. 


(38b) 


We  can  also  develop  another  form  from  (37)  for  w  »  Aw  by  retaining  the 
exp(-x/Km2L02)  but  setting  (l+x)“^  «  x"*!.  Equation  (37)  is  then  reduced  by 
utilizing  the  Laplace  transform  (4.5.29)  of  Reference  [9J,  and  by  introducing 
the  microscale-governed  frequency  AD2  -  16Km2AU2/3.  We  obtain  (for  q  >  1/2), 


ryj—  ( k  L  w\^2  q 

V“-»  ■  I ?  \rsj  Vi/2(2"/ia) 

'  S  (sr)  f°r  “  «  at!  «<>■=> 

'  <W  L  (ft)  q  «P<-2»/AO)  for  t>  »  40. 
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Here,  Kq_i/2(2w/Aft)  is  a  modified  Bessel  function  in  the  notation  of  Reference 
[10],  and  the  two  asymptotic  forms  are  easily  obtained.  We  note  that  the 
second  form  in  (38c)  agrees  with  the  second  form  in  (38b)  thus  confirming  the 
validity  of  the  approximations.  It  is  also  clear  that  the  power  spectrum  cuts 
off  abruptly  as  u>  approaches  and  exceeds  Aft  ^  AU/iQ. 

Just  as  in  case  B,  we  note  that  the  power  spectrum  starts  off  flat  below 
w  »  AU/L0,  and  then  it  decays  as  w"2/3  until  the  cut-off  at  «  *  AU/io. 

D.  SUMMARIZED  RESULTS  AND  APPROXIMATED  POWER  SPECTRUM  FOR  GENERAL  CASE 

The  results  found  in  this  section  amount  to  the  following.  We  assume  a 
random  velocity  distribution  with  mean  UQ  and  variance  AU?.  The  component  of 
mean  velocity  across  the  beam  is  U?>  We  have  defined  a  power  spectrum  of 
focal-spot  radius  fluctuations  that  is  proportional  to  the  power  spectrum  of 
angular  fluctuations.  The  latter  is  named  Wg(u>).  It  has  to  be  a  function  of  Uj 
and  AU^.  We  have  defined  some  frequencies  by  dividing  these  velocities  by 
length  scales: 


°T 


VLo 

,  Aw  *  16AUZ/3L 

o 

(39) 

k«UT 

,  Aft2  =  16k2AU2/3  . 
in 

Furthermore,  we  have  found  that  Wg(w)  is  the  difference  of  two  parts,  see 
Equation  (33).  In  Tables  I  and  II  we  summarize  the  results  for  Wg(u>,q), 
bearing  in  mind  that  we  subtract  the  q  -  11/6  equations  from  the  q  -  5/6  ones 
to  get  the  power  spectrum. 


It  is  interesting  to  compare  the  random-flow  with  the  frozen-flow  results. 
The  spectra  are  very  similar.  We  would  like  to  suggest  an  approximate  formula 
for  engineering  purposes.  To  do  so,  we  note  from  (31)  two  properties: 


V°>  '  7  / 


d“V“>  -  Tbr  ir r<11/6)  W1'3' 

o 


(40) 


■[ 


32R6(t) 


3t 


T«0 


w 

1  C.  2..  ,  ,  15.7  =  L  r(7/6)  ,  .  ,7/3  .  2  .  .  2, 

7  J  d““  »,<«>  -  isr  •  i — y^-  <*my  ("t  +4“  > 


The  second  property  expresses  the  variance  width  of  the  power  spectrum.  It  is 
equal  to  the  mean  square  width  because  Wg(u)  is  symmetric  around  w  ■  0  so  that 
we  utilize  cosine  transforms  in  the  preceding.  The  first  property  simply 
recalculates  the  angular-fluctuation  variance  and  states  that  it  is  the  zeroth 
moment  of  Wg(w).  The  interesting  point  is  that  (40)  is  a  function  of  uj2+Auj2. 

We  already  know  that  the  flat  part  of  Wg(u)  (for  u  <  avp,  Aoj)  does  not  contribute 
appreciably  to  Rg(0)  because  it  can  be  seen  readily  from  (31)  if  we  set  1+x  ■  x 
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TABLE  I 

Frozen  Flow  (AU=0) 


I 


I 


that  the  resulting  power  spectrum  then  has  no  flat  part,  and  we  have  already  seen 
that  this  is  a  valid  approximation  in  computing  Re(0).  Consequently,  the  flat 
part  is  even  less  important  in  computing  the  second  of  (40).  Since  most  of 
the  area  under  W q(u))  lies  in  the  m“2/3  portion  (q  ■  5/6  yields  an  to- 2/3 
behavior  that  dominates  over  the  q  ■  11/6,  m“8/3  part),  it  appears  to  us  that 
this  portion  of  the  curve  is  governed  by  a  parameter  m*/ (mpj2+Am2)  .  Consequently, 
we  suggest  that  the  following  formula  is  useful  over  the  entire  range  of  m: 


W0(«)  = 

U0(u,it) 


15.7  c2L 
8tt  L 

o 


[W0(a»,5/6)  -  WQ(m,ll/6)], 


1/2-q 

exp[-m2/(fi2+A02)] , 


x  U(l/2,  3/2-q,  (u)2+u)T2+Aa)2)/(n2+Afi2)]  . 


(41a) 


This  formula  clearly  satisfies  (40)  because  it  reduces  to  (36a)  if  we  set 
Am  -  0,  and  (36a)  yields  (40)  with  Am  =  0.  It  is  obvious  that  by  substituting 
w<j»2+Au)2  for  in  (36a)  we  satisfy  (40)  in  general,  and  we  obtain  the  above 
formula  (41).  When  oj2  «  +Afl2f  we  can  utilize  the  approximation 


W0(u),q)  * 


r(q-l/2) 

r(q) 


1/2 


(41b) 


Comparing  (41b)  to  the  results  in  Tables  1  and  II,  we  note  that  we  obtain  the 
previous  results  to  reasonable  approximation.  If  we  allow  m-p  •+■  0,  we  note 
that  the  (41b)  approximation  is  in  error  for  m  «  Am  by  a  factor  Jn,  and  by  a 
factor  T  (5/6)  1.1  for  m  »  Am.  This  appears  sufficient  for  engineering 

purposes.  Furthermore,  (41b)  is  in  error  when  ivp  -*■  0  for  m  %  Afi,  but  here  the 
spectrum  is  very  small  and  its  precise  shape  is  governed  by  the  dissipation 
range  of  turbulence  for  which  the  precise  shape  is  not  crucial  anyway.  By 
choosing  (41)  as  we  have,  we  at  least  assure  ourselves  of  the  proper  normaliza¬ 
tions  (40)  although  the  cut-off  shape  at  m  ^  Afi  is  not  described  properly.  The 
transition  regions  are  probably  also  imprecise. 

Therefore,  in  conclusion  it  appears  to  us  that  the  power  spectrum  of 
angular  fluctuations,  hence  that  of  focal-spot  scintillation,  is  adequately 
described  by  (41) .  That  is  to  say, 
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2  2 

(i)  the  spectrum  is  flat  for  w  <  (ui^,  +Aw  ). 

2  22  “1/3 

(ii)  it  falls  off  as  [ui  /(wT  +Aw  )] 

for  (w  2+Aw2)1/2  <o i<  (fl  2+Afl2)1/2. 

T  1  2  2  1/2 

(iii)  it  cuts  off  sharply  beyond  w  *  (0T  +AJ2  ) 


To  our  knowledge,  no  calculations  of  We(w)  exist  in  the  literature  for 
this  general  case.  However,  Lawrence  and  Strohbehn[ll]  cite  untranslated  work 
by  A.  S.  Gurvich,  M.  A.  Kallistratova,  and  N.  S.  Time  (Izv.  Vyssh.  Ucheb. 
Zaved.  Radiofizika  11,  pp.  1360-1370,  1968)  for  the  Taylor's  hypothesis  case 
(Aw  -  0).  Figure  10  of  Reference  [11]  shows  a  plot  of  essentially  w  Wq(w)  vs. 
w  which  appears  to  be  misinterpreted  in  the  text  of  Reference  [11].  All  three 
portions  (i)  -  (iii)  of  the  behavior  mentioned  above  are  borne  out  by  this 
graph  which  appears  to  us  to  correspond  to  wx  ^  1  Hz  (note  that  WT  "  2ir  f T) , 
and  fix  *  300  Hz  (because  the  peak  of  the  curve  is  at  icmL0fx  «*  50  Hz) . 


APPENDIX  I 


THE  VARIANCE-WIDTH  THEOREM  FOR  FOCAL-SPOT  AREA 


The  purpose  of  this  Appendix  is  to  derive  (3)  in  a  very  general  fashion 
and  to  stress  the  approximations  in  it.  These  do  not  appear  to  be  understood 
as  clearly  as  is  desirable  in  context  with  the  derivation  in  Reference[l] . 

Referring  for  notation  to  Reference[l] ,  and  specifically  to  Section  III, 
we  note  that  the  irradiance  in  the  focal  plane  (for  a  focussed  laser  beam)  is, 

I(r)  «  |  BOr,^)  U^)  exp(ik  p -^/L)  |  2  (A-l) 


Here,  ^(r^  is  the  transmittance  amplitude,  the  so-called  pupil  function,  of 
the  aperture  at  z  =  0,  and  coordinate  r^  =  (p^O).  Previously,  UQ  was  chosen  to 
be  Gaussian.  At  present,  it  can  be  considered  any  generally  accepted  shape. 

The  "sharp-beam"  choice  UQ(ri)  -  1  for  px  i  r„,  and  zero  otherwise,  yields  the 
Airy  disk  in  (A-l)  when  there  are  no  atmospheric  effects  (B  =  1) . 

The  approximations  in  (A-l)  are: 

(i)  radiation-zone,  kL  »  1 

(ii)  near  field  of  the  aperture,  L  «  kr  2 

(iii)  sagittal,  krQ4  «  L-*  ° 

These  are  not  very  great  restrictions  at  all  for  X  =  10.6  ym  and  higher- 
frequency  propagation.  At  10.6  ym,  one  finds  kr02  ^  100  km,  and  (kr04)l/3  30  m. 

Clearly  1  km  <  L  <  10  km  yields  correction  terms  to  (A-l)  that  are  quite  small. 

To  proceed  further,  we  note  that  the  expression  inside  modulus  bars  in 
(A-l)  is  a  two-dimensional  Fourier  transform  of  a  function  of  p1.  The  conjugate 
variable  is  kp/L  »  q,  and  consequently  the  irradiance  is  a  function  of  L  and 
of  q.  Therefore 


!(q)  a  I  J'd2P1  B(r,r1)  Uq (rx)  exp(iq*p1)  |2.  (A-2) 


When  B(r,r^)  =  1,  we  obtain  Ic(q),  the  free-space  irradiance  pattern  in  the 
focal  plane  z  =  L. 

Only  one  more  approximation  will  be  made:  it  is  assumed  that  every  vector 

“►  “►  -  I  -*■  I  < 

r~ri  tor  |r-r^  |  rQ  is  perpendicular  to  the  plane  z  =  0.  This  yields, 

(iv)  small-angle  approximation,  r  <<  L 
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This  approximation  yields  extremely  small  errors  and  is  easily  satisfied  in  all 
practical  situations.  It  is  needed  only  so  that  <B(r ,ri )B*^?,?2) >  can  be  re“ 
placed  by  the  mutual-coherence  function  (mcf)  with  axis  r  _j_  r r 2  • 

Because  <  B(?,r^)B*(r ,1*2)  >  is  a  function  of  the  difference  coordinate 
tp  -  r^-^2*  1  Fourier-transform  pattern  IB(<J)  can  be  defined  as  follows 

<IB(q)  =  Ji d2Ap<B(r,r1)  B*(r,r2»  exp(iq*Ap).  (A-3) 

The  average  of  Equation  (A-2)  can  be  rewritten  as  a  convolution  integral 
of  lg  and  IQ  because  of  this  property  of  the  mcf  of  the  normalized  spherical- 
wave  field  B: 

<IB(^)>  ““2  /d2q  1  VV  (A-4) 


This  formula  for  the  average  irradiance  is  not  new,  but  it  is  useful  to  stress 
its  convolution  property.  It  states  that  the  free-space  irradiance  pattern 
lQ(q)  is  broadened  by  the  pattern  Ig(q)  which  depends  only  upon  the  mcf  of 
normalized  spherical-wave  fields  in  turbulent  air. 

2  2 

We  now  define  variance  widths  for  <I>,  IQ,  and  IB,  namely  a  ,  aQ  ,  and 
a  2  respectively: 

D 


/d2q  q2  <  I (q)  > 

Jd2q  q  <  I  (q)  > 

/d2q  <  I (q) > 

Jd2q  <  I (q)  >  _ 

(A-5) 


2  2 

Likewise,  definitions  for  oQ  and  oB  can  be  written  down,  although  omitted 
here.  If  (A-4)  is  inserted  into  numerators  and  denominators  of  (A-5),  one  finds 
the  variance-width  theorem  after  some  rearrangement  of  terms: 


2 

a 


(A-6) 


Note  that  the  derivation  depends  only  upon  the  validity  of  (A-2),  which 
requires  assumptions  (i)-(iii)  and  the  fact  that  the  mcf  <  B(1)B*(2)  >  (in 
shorthand  notation)  is  a  function  of  r^-?2*  applying  (A-6)  to  effective 
irrad lance-pat tern  widths,  it  is  also  noted  that  each  o2  need  be  divided  by 
(k/L)2  to  yield  actual  widths.  Thus,  if  we  define  ry  ■  (L/k)o,  r.  ■  (L/k)o  , 
and  r. _  ■  (L/k)o_,  we  obtain  °  0 

L  D  D 
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(A-7) 


Lo 


+  r. 


LB 


respect  to  the 
for  the  first 
,  rLn2  are 

defined  likewise  by  replacing  <  I$,L) >  by  IQ(£,L)  and  Ig(p,L)  ■  Ig(q)  respec¬ 
tively.  Finally,  It  Is  noted  that  an  expression  for  r. _2  is  readily  obtained 
from  the  definition, 

rLfi2  -  /d2p  p2  IB(q)  / / d2p  IB(q),  (A-8) 


For  cylindrically  symmetric  beams  ql  (q)  is  antl-symmetric  with 
origin,  hence  the  first  moment  of  I  ¥<f)  is  zero.  The  same  holds 
moment  of  In(q).  Hence  rt  is  prooerlv  defined  by  (1).  and  rT „ 


by  utilizing  Equation  (10)  of  Reference[l]  for  the  mcf.  Unfortunately,  part  of 
the  exponent  in  that  equation  is  in  error.  The  correct  form  for  <  B(1)B*(2)  > 
with  $  -  ri~?2  is. 


< B(1)B*(2)>  -  exp 


(A-9) 


When  this  is  inserted  into  the  definition  of  lg(q) ,  which  is  then  inserted  in 
turn  into  (A-8) ,  one  obtains  via  a  development  analogous  to  that  in  Section  IV 
of  Reference [1 ] , 


T  3  2 
L  e 

24tt 


oo 

/ 


dKK3$(K) 


1.3  L3k 


1/3 


m 


(A-10) 
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APPENDIX  II 


THE  RANDOM-VELOCITY  HYPOTHESIS 


The  difficulty  in  (29)  is  that  the  autocovariance  <  fie  (r^t)  fie  (r2;t  +  t)> 
depends  on  the  time  separation  t.  We  will  discuss  a  simple  hypothesis  that 
enables  us  to  relate  this  autocovariance  to  the  t  ■  0  case.  Consider  5e(r;t) 
in  turbulent  air,  it  is  proportional  to  the  air-density  fluctuation  6N(r,t). 

The  continuity  equation  states  that 


dN/dt  +  NV#U  -  0 


(A-ll) 


where  N  =  N0  +  5N,  and  5  is  the  gas  velocity.  Thus  (A-ll)  describes  the  time 
change  of  fiN  because  we^may  assume  that  dN<j/dt  *  0.  In  a  hydrodynamical^d in¬ 
scription,  both  5N  and  U  are  functions  of  r  and  t,  and  dN/dt  *  3N/  3t  +  U*vN. 

The  second  term  of  (A-ll)  is  related  to  expansion  or  compression  of  the  gas  and 
the  first  term  simply  expresses  conservation  of  density  in  a  unit  parcel  of  air 
as  it  flows  from  one  location  to  another  (perhaps  changing  its  shape) . 

If  we  assume  that  "sources"  of  turbulence  <Jo  not  lie  directly  in  or  close 
to  the  propagation  path,  then  we  may  assume  Nv»U  «  0  and  consider  the  turbulence 
adequately  described  by  dN/dt  -  0,  i.e.,  by  flow  governed  by  a  wind-velocity 
vector  ft(r,t).  It  is  probably  not  a  very  wrong  assumption  to  thus  ignore  com¬ 
pression  and  expansion  of  the  gas,  and  it  simplifies  matters  because  it  implies 
that  6N(r,t)  can  be  written  in  general  as 

t  +  T 

6N(r,t)  -  fiN (r  +  J dt'  U(r,t  ’),t)  (A-12) 

t 

and  similarly  for  fie  «6N.  It  is  easily  seen  that  (A-12)  implies  dN/dt  -  0  by 
allowing  x  0  and  keeping  the  first-order  Taylor-series  terms.  These  yield 
dfiN/dt  +  U*^6N  -  0  which  is  identical  to  dN/dt  *  0  because  Nq  does  not  vary. 

We  can  apply  (A-12)  to  (29)  to  obtain: 

t  +  x 

fie(r2;  t  +  x)  -  fie(r2  -  j  dt'  U(r2,t');  t)  (A-13) 

t 

If  this  is  inserted  into  (29)  we  obtain  after  a  coordinate  transformation: 

t  +  x 

jjLK'Ar '  +  J  dt'  U(r2,t')J  (A-14) 

4>^(ifc;x)  -  <y*d^Ar'fie(r1' ;t)  fie2';t)  e  t  > 
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f 


where  Ar'  ■  ri'-^',  rj'  ■  r^,  and  ^  is  the  spatial  coordinate  in  the  right- 
hand  side  of  (A-13) .  The  coordinate  A?'  is  time-dependent,  but  since  that  time 
dependence  shows  up  only  in  the  boundaries  of  integration,  effectively  replaced 
by  +  00  in  all  dimensions,  we  can  disregard  it. 


The  cardinal  assumption  to  be  made  at  this  point  >is  that  the  velocity  ; 
structure  is  statistically  independent  of  the  spatial  structure  of  6e. 

It  allows  us  to  write  t  +  t  ' 


(K;x)  - / d3Ar'  <  6e(r1')6e*(r2')>  e 


iK*/£r *  lit*  f  dt’U(r2,t') 

(■  '■  '  ) 


, (A-15) 


The  space  dependence  of  U  has  not  been  written  down  in  the  last  form  because  an 
ensemble  averag£  over  space  and  time  is  implied  by  the  brackets.  It  is  also 
plausible  that  U(t')  does  not  vary  locally  very  much  with  t'  in  a  time  interval 
t  £  1  sec.,  or  t  less  than  soijie  other  short  time  in  which  4>4(K;t)  Is  essentially 
zero.  Furthermore,  velocity  U  is  most  likely  to  consist  of  a  steady  wind  com¬ 
ponent  and  a  random  component,  l.e,,  its  distribution  is  probably  adequately 
described  by  a  Gaussian  with  non-zero  mean.  These  two  statements  imply,  1 


<exp 


iiE-  f 

“  «- 


t  +  T 


dt'  U(t') 


iK-U  x  2 

e  0  exp (-2*  Ain  x  ) . 


1  (A-16) 


Here,  UQ  is  the  mean  compnent,  and  we  note  that  only  its  component  U<j.  in  a 
plane  normal  to  the  propagation  direction  occurs  in  the  first  exponential. 
Likewise,  AU>[  is  the  variance  of  wind  velocity  in  this  plane.  If  the  Wind 
variation  is  isotropic,  then  it  follows  that  AU^  ■  (2/3)AU^.  Finally,  all  this1 
is  summarized  by  inserting  it  into  (A-16)  which  is  inserted  in'  turn  into  (A-15) 
to  obtain 


$4(K;t) 


1K-Utx 


4  2  2  2 

exp(-  y  K  AU  x  )  4>(K) 


1 

(A-17) 


This  set  of  assumptions  leading  to  (A-17)  is  also  known  as  Favre's  hypothesis [8 ] . 
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